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ABSTRACT 

We examine the stability and observational consequences of mixing induced 
by ^He burning in the envelopes of first ascent red giants. We demonstrate that 
there are two unstable modes: a rapid, nearly adiabatic mode that we cannot 
identify with an underlying physical mechanism, and a slow, nearly radiative 
mode that can be identified with thermohaline convection. We present observa- 
tional constraints that make the operation of the rapid mode unlikely to occur 
in real stars. Thermohaline convection turns out to be fast enough only if fluid 
elements have finger-like structures with a length to diameter ratio l/d > 10. 
We identify some potentially serious obstacles for thermohaline convection as the 
predominant mixing mechanism for giants. We show that rotation-induced hori- 
zontal turbulent diffusion may suppress the ^He-driven thermohaline convection. 
Another potentially serious problem for it is to explain observational evidence of 
enhanced extra mixing. The ^He exhaustion in stars approaching the red giant 
branch (RGB) tip should make the ^He mixing inefficient on the asymptotic gi- 
ant branch (AGB). In spite of this, there are observational data indicating the 
presence of extra mixing in low-mass AGB stars similar to that operating on the 
RGB. Overmixing may also occur in carbon-enhanced metal-poor stars. 

Subject headings: stars: abundances — stars: evolution — stars: interiors 



1. Introduction 



There is strong observational evidence for deep mixing in the radiative envelopes of 
low-mass (M < 2 Mq) red giant branch (hereafter, LM-RGB) stars. Changes in light ele- 
ment abundances (such as Li, C, N) and in the ^^C/^^C ratio as a function of luminosity have 
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) . The observed pattern requires at least a component of in situ mixing. This 



extra mixing could have consequences for other species (such as ^He) that are not directly 
observed. Indeed, in spite of the predicted efficient production of ^He in low-mass main 
sequence (MS) stars, its G alactic abundance has been nearly constant since the epoch o f Big 



Bang nucleosynthesis (e.g. , lTosilll998l : iBania. Rood, fc Balserll2002l : IVangioni et al.ll2003l ) . To 



explain this, semi-empirical stellar evolution models have shown that the carbon depletion 
due to extra mixing in LM-RGB stars should unavoidably be accompanied by a strong ^He 



destruction t hat counterbalances its production on the MS (IRood, Bania. fc Wilson 



Hoganl I1995I : ICharbonnell Il995l : IWeiss. Wagenhuber. fc Denissenkovl Il996l ) . 



1984 



However, the true physical process that is responsible for mixing has resisted identifica- 
tion. Rotation ally induced mixing has be en an implied underlying mechanism since the pio- 
neering work of ISweigart fc Mengell (119791 ). but there a re serious difficulties in reconciling the 
observed mixing pattern w ith theoretical predictions (IChaname. Pinsonneault. &: Terndrup 
2OO5I : iPalacios et al.ll2006l ). Moreover, it can be shown that rotational mixing solely de- 
pendent on the local angular velocity gradient q = (dlnfl/dlnr) (e.g., shear mixing) will 

be self-quenching. Indeed, the empirically constrained mixing rate fmix ^ 10~^-10 

4 , 



-2 —1 

cms 



(see §[5]) is faster (as it should be) than the mass inflow rate |f| < lO^^cms^^ in the radiative 
zones of LM-RGB stars. Hence, if fmix were proportional to {n > 0) then the flattening 
of the rotation proflle by the accompanying angular momentum redistribution would quench 
the mixing very quickly. 



A very different class of solutio n has recently been proposed bvlEggleton. Dearborn, fc Lattanzio 



(120061 ) (hereafter, EDL06; see also lEggleton. Dearborn, fc Lattanzio! 120071 ). While investi- 
gating the core He flash in a low-mass ni odel star near the RGB tip with the code D j ehut 
Dearborn. Lattanzio. fc EggletonI (12006! ) noticed some gas motions in the radiative zone 
above the H-burning shell, in addition to convective motions driven by ^He burning in 
the core. In their follow-up papers, they have made a conclusion that these additional gas 
motions are due to the Rayleigh- Taylor instability (RTI) driven by ^He burning. They have 
noted that, unlike other nuclear reactions in stars, the reaction ^He (^He, 2p)^He decreases the 
mean molecular weight fi. EDL06 have argued that, even though the decrease of /i is minute 
(|A/i| ~ 10"'^), it is the resulting local inversion of the /i-gradient = {din fi/dln P) < 
that has led to the RTI in their 3D simulations. They claim that velocities of gas motions in- 
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duced by the RTI are "comparable to the velocity of the normal convection" [vc ~ 10^ cms~^) 
and that this is consistent with the simple heuristic estimate of ~ gHp{Afi/ fi), where g 
and Hp are the local gravity and pressure scale height. 

A yU-gradient inversion from ^He burning is a predicted consequence of standard stellar 
evolution. After a low-mass star has left the MS its convective envelope first deepens (the 
first dredge-up) and then its bottom starts to recede. At the depth of its maximum inward 
penetration the convective envelope imprints a discontinuity in the chemical composition 
profile. Later on, the advancing in mass H-burning shell will erase this discontinuity. During 
this event the evolution of the star slows down, which produces bumps in differential luminos- 
ity functions of stellar clusters, and the star itself makes a tiny zigzag on the HRD diagram 
(Fig. [H upper panel). However, before the major H-burning shell erases the composition 
discontinuity, a shell in which ^He burns down and which advances in front of the H-burning 
shell will cross the discontinuity first. The ^He(^He,2p)^He reaction reduces the mean molec- 
ular weight locally. However, this reduction is so minute, A/i ^ ^"^AX^/Q ~ — 10~^ (here, 
fi ~ 0.6, and is the ^He mass fraction that can reach up to a value of ~2 x 10^'^ in the 
envelope of a low-mass star ascending the RGB), that it is not seen on the /x-profile until 
the II depression will find itself in the chemically homogeneous region of the radiative zone 
previously occupied by the convective envelope. This happens close to the bump luminosity 
(Fig. [I]). 

In this paper we investigate the ^He instability in more detail. First of all, we note 
that the RTI would not be expected in compressible, and hence stratified, stellar material. 
However, a /x-gradient inversion may trigger instabilities related to convection. In section 
2, we discuss general criteria for convective instability taking into consideration that, over 
longer timescales, fluid elements can exchange heat with their surroundings, thus modifying 
the background temperature stratification. We identify two families of solutions: a rapid 
mode with a nearly adiabatic thermal structure and a slow mode with a nearly radiative 
structure. The EDL06 results appear to correspor id to the rapid mode, while the slow mode 



can be identified as thermohaline convection (e.g.. I Vauclairll2004j ) . On the basis of previously 
published diffusion coefficient estimates it is likely to be triggered in the envelopes of red 
giants. However, there are significant (and uncertain) assumptions related to the actual 
operation of this instability. We discuss key features determining the efficiency of the slow 
mode: the assumed geometry of the fluid elements, which directly impacts the timescale for 
exchanging heat; the potential impact of horizontal turbulence in suppressing the instability; 
and the predicted depth of mixing. In section 3 we demonstrate that the alternative approach 
based on the linear stability analysis of the underlying conservation equations also leads to 
the conclusion that thermohaline convection may be suppressed by the horizontal turbulent 
diffusion. 
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In section 4 we evaluate the impact of any proposed mixing mechanism on the thermal 
structure of the red giant branch stars. We argue that the rapid mode, which is similar 
in its properties to the originally published results of EDL06, can be ruled out because it 
would induce strong feedback on the thermal structure of giants and would predict a mixing 
pattern contrary to observations. We also demonstrate that both the slow mode and prior 
empirical mixing estimates would not disturb the thermal structure of giants. In section 4 
we assess the overall promise of '^He mixing as a mechanism. We find that it may be an 
attractive solution, but identify several potentially serious drawbacks on both observational 
and theoretical grounds. In particular, we argue that previously published estimates of 
horizontal turbulence would be sufficient to suppress the instability and that the naturally 
expected depth of mixing and trends with luminosity may be in conflict with observational 
trends. 



2. General Criteria for Convective Instability 

In the presence of a positive /x-gradient dfi/dr (a negative V^), a fluid element displaced 
vertically upwards will find itself surrounded by material with a higher /i. Whether the fluid 
element will continue to rise depends on how efficiently it can exchange mass and heat with 
its surroundings. It is easy to anticipate that heat diffusion will favor the instability by 
reducing the temperature gradient inside the rising fluid element V = (ci In T'/ci In P), as 
compared to its adiabatic value Vad, and shghtly increasing the gradient in the surrounding 
medium V = (dlnT/dlnP), as compared to its radiative value Vrad in the absence of mixing. 
In a radiative zone, we always have Vrad < V < V < Vad- On the other hand, molecular 
diffusion and horizontal turbulent diffusion (if the latter is present it will also contribute to 
heat diffusion) will decrease the //-contrast between the fluid element and its surroundings 
|A/x| = |/x'(r) — /i(r)|, thus hindering the development of convection. Current theories of 



rotat i onal mixing predict the existence of strong ho rizontal turbulence (IChaboyer fc Zahn 



19921 : iMaederl l2003l : iMathis. Palacios. &: ZahnI 120041 ): even if the vertical turbulence is not 



strong enough to produce mixing such a horizontal turbulence could strongly impact an 
instability driven by differences in composition. Let us carry out a general investigation of 
effects produced by these diffusion processes on the convective instability of the radiative 
zone in the presence of a /i-gradient inversion. The magnitude of these effects depends on 
the Peclet number that represents a ratio of thermal and mixing time scale 

„ Tth Iv f d^ ^ 

Pe = ■ 



(K + D^) \l 
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where rth ~ (P/{K + D^a) and Tmix ^ l/v. Here, we consider a possibility that the convective 
motions may be organized in elongated narrow structures resembling "fingers" whose length 
/ is much larger than their diameter d. In this case, the fiuid element velocity v, thermal 
diffusivity K, and horizontal turbulent diffusivity should appropriately be averaged over 
the finger length scale /. For simplicity, we will represent a "finger" by a spherical fiuid 
element of the diameter d that travels the path I > d before it dissolves. Modeling a "finger" 
with a cylinder would introduce factors of order unity in our derived relations. In the limit of 
a high Peclet number, mixing is so fast compared to heat exchange that T' undergoes nearly 
adiabatic changes, hence V ~ Vad- In this limit, when the rising fiuid element dissolves 
it has a lower temperature than its surroundings therefore it reduces T locally making V 
steeper and closer to Vad- In the opposite limit of a low Peclet number, mixing is so slow 
that the radiative and turbulent heat fiux from the surroundings can warm up the fiuid 
element thoroughly as it rises. This brings V close to Vrad while V remains nearly equal to 
Vrad because the fiuid element will have T' ^ T when it dissolves. 

In our further anal ysis we will us e sim ple heuristic relation s obtained in the mixing 



length approximation by iMaederl (1l995l ) and iTalon fc ZahnI (1l997f ) in their investigations of 
shear instabilities in rotating stars. We admit that the radiative zone may host both the '^He- 
driven convection and some other extra mixing of a nonconvective origin (e.g ., rota t ion-dr iven 



turbulent diffusion or meridional circulation) at the same time. Following IZahnI (Il992l ). we 
assume that rotation-driven turbulence is highly anisotropic, with horizontal components of 
the turbulent viscosity strongly dominating over those in the vertical direction, ^ D^. 
We consider a process with a diffusion coefficient -Dmix = vl/3. At the present stage we leave 
-Dniix unspecified, and simply solve for the coupled diffusion of heat and chemicals to evaluate 
the range of diffusion coefficients over which an instability occurs. In the next section we 
compare with specific (and previously published) estimates of diffusion coefficients. 

The degree to which convection modifies the thermal structure of the radiative zone 
depends on its heat transport efficiency 



AEtr.., V'-V A.ix Pe 



AEex Vad - V 2{K + Z^h) V V 6 



'1^ 



The quantity F measures the deficit or excess of energy transported by rising or sinking 
turbulent elements Aii^trans with respect to the energy Aii^ex the elements gain or lose through 
the radiative (K) plus turbulent (Dh) heat exchange with the surroundings. 

For an ideal gas, the square of the Brunt- Vaisala (buoyancy) frequency is 

iV2 = -^(V'-V + V,-V;). (2) 
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This equation takes into account that horizontal diffusion (molecular plus turbulent) may 
change the mean molecular weight of the fluid element yu' during its motion, so that 7^ 0. 
The convectiv e instability se t s in w hen N"^ becomes negative. By analogy with equation ([T]) 
and following iTalon fc ZahnI (119971 ). we introduce a /i-transport efficiency 



3D. 



v;, 2(z/^oi + Dh) 



(3) 



where z/moi is the molecular diffusivity. Supplementing the radiat ion luminosity with a con- 
tribution to heat transport by convection in the same manner as iMaederl ( 1l995l ) dealt with 
shear mixing, we obtain the following relation between V, Vrad and Vad^ 



Vrad + 6 



r+1 K 



ad 



1 + 6^-^+^'- flV 



(4) 



Combining eqs. (P-Hl), the instability condition A^^ < can be transformed into 



- V^ > 



i + r + eP^Mh (1)' 



(Vad — Vrad) 



(5) 



2.1. The Adiabatic Limit 

In the case of F ^ 1, equations ([1]) and (jl]) give V' ~ V ~ Vad, which we call "the 
adiabatic limit" . From equations ([1]) and ([3]) it follows that 

r. = 3;i^r, (6) 

hence F^ 1 as soon as F ^ 1 because i/moi ^ K (Fig. [2]). Using these constraints, the 
condition ([5]) is transformed to 

D^^ >\k ^^^r=^ » K. (7) 

Radial displacements / ~ 10^ cm of fluid elements with velocities of order 5 x 10^ cms~^ 
observed by EDL06 in their 3D red giant model above the major H-burning shell correspond 
to a diffusion coefficient -Dmix ~ \vl ~ 2 x 10^^ cm^ s~^ ^ K. The same or even higher order 
of magnitude estimate can be obtained if we calculate -Dmix ~ vHp, where ~ gHpAfi/ fi, 
as proposed by EDL06 (see next section). So, it appears that the ^He-driven mixing in the 
3D red giant model of EDL06 somehow wound up in the metastable adiabatic limit. We 
explore the consequences of such a rapid mixing process for both surface abundances and 
the thermal structure in section 4. 



2.2. The Radiative Limit 



Let us now consider a situation wlien T{l/dY <^ 1- Tliis also means tliat F ^ 1 because 
we have assumed that d < I. In this case, equations ([1]) and (jl]) give V' ~ V ~ Vrad, 
therefore we will refer to it as "the radiative limit". Given that in equation ([6]) the ratio 
K/vrnoi ^ 1 (solid curve in Fig. [2]), the assumption that F <^ 1 does not automatically lead 
to F^ < 1 unless > K. Hence, we have to consider the cases of F^ 1 and F^ < 1 
separately. 

a). In the radiative limit (F <^ 1), values of F^ ^ 1 can be met only if D]^ -C K. Using 
these constraints, the condition ([5]) is simplified to 

Ml 



^mix < ^ ^ IV 

Vad — Vrad 



d 



The right-hand side of ([8]) a dequately reproduces both the diffusion coeffi cient for thermo- 
haline convection derived by iKippenhahn. Ruschenplatt. fc Thomas! (jl980l ) 

Z^Kipp = 3-^^— IV^I, (9) 

Vad — Vrad 

who argued that / sh ould be of ord er d, and the rate of mixing by elongated narrow "fingers" 
(/ > d) advocated by luirichi Jl972h 

Dvmc^=\'k \ ^ . (10) 

3 Vad -Vrad V"/ 

Thus, in the radiative limit with <^ K we can readily identify mixing with thermohaline 
convection. Substituting expressions (P-[TOl) in place of i^mix, we find that, for thermohaline 
convection driven by the ^He burning, 

r — f-Y -Pmix / dV |V^| 

'"^GKyiJ 2K [ij Vad -Vrad ' ^ ' 

as we assumed. 

b). If is not negligibly small compared with K then we sure have ^ t'moh ^ind 
relation can be re-written as 

r. = 3^r. (12) 



For F^ < 1, the condition takes a form 
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(compare it with condition 5 from IVauclairl |2004| ). In the radiative zone of an LM-RGB 
star, the left-hand side of f|T3|) is of order 10~^ at most. A profile of the quantity K in 
the radiative zone of our 0.83 bump luminosity model is plotted with dashed curve in 
Fig. [21 Given that K < 10^ cm^ s~^ for r < 0.1 i?© and comparing the ratio Dy^/ K with 
the number 10~^, we conclude that the horizontal turbulent diffusion with > 10^ cm^ s~^ 
may hinder the dev elopment of convectiy e instability. Interestingly that such values of 
have been found by lPalacios et al.l (120061 ) in their low-metallicity 0.85 bump luminosity 
model even for the less favorable case of solid-body rotation of the convective envelope. For 
the case of differential rotation of the convective envelope, which is suggested by observed 
fast rotation of horizontal branch stars, they have obtained D^i > 10'' cm^s~^. For such 
large values of Dh, the right-hand side of flT3l) exceeds the expression on the left-hand side at 
least for r < 0.1 i?0 (Fig. [JJ, therefore the convection (thermohaline mixing) will probably 
be suppressed there. 



Ulrich's Solution in the Presence of Horizontal Turbulent Diffusion 



In this section we demonstrate that lUlrichl (119721 ) could have come to a conclusion 
about the suppression of thermohaline convection by the strong horizontal turbulent diffusion 
similar to that made by us (condition [T^ if he had included in his equations. To do 
this, we start with the linearized equations of conservation of momentum, thermal energy, 
and chemical composition similar to those used by Ulrich (his equations 1-3) but with the 
diffusion coefficient taken into account 



dv 
2 — 
dt 

d6\nT 

dr~ 

d 6 In fi 
dt 



V 



-:V — (yf^lnT+fyf^lnU, 

d^ 



rf2 



51nT, 



E, 



V b In /i. 



d2 



(14) 
(15) 
(16) 



Here, v = z/moi+^'rad is the total (molecular plus radiative) viscosity, SliaT = In T'(r)— In T(r), 
Slnfi = Inyu'(r) — In /i(r), other symbols having been defined previously. 

The characteristic polynomial for the linear system of ODEs (lT^ - lT6ll can be written in 
the following form: 

{cun^f + [1 + (rth/r,) + (nh/r^)] (cuTth)' + (17) 
[(rth/r.) + (rth/r,) (nt/r^) + (nh/r^) + {rlN^)] (ujnt) + 
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[(rth/r.) (rth/r^) + (r^liV^ + {njr,) (r.liV^)] = 0, 

where oj = ^nr^l, is the eigen frequency of stable {u < 0) or unstable {uj > 0) mode, 
= 5'(Vad— Viad) Hp^ and A*"^ = gVfj, Hp^ are the squares of the T- and /i-component of the 
Brunt-Vaisala frequency for the ideal gas, N'^ = Np+N"^, while = d'^/i', Tth = d"^/ {K+Dy,), 
and = (i^/(t'moi + D^) denote the viscous, thermal, and horizontal diffusion timescales for 
the fluid element. Our equation (fT7|) is equivalent to Ulrich's dispersion relation (10). 

The only real root of the polynomial (IT7I) is plotted as a function of (D^/K) in Fig. [3] 
for the ratio d/Hp = 0.01 and three different values of A*"^ = —10^^ (solid line), — 10~^ 
(short-dashed line), and — 10~® (long-dashed line). The first value of A^^ is close to the 
minimum one found in the region of the /i-gradient inversion produced by '^He burning in 
a low-metallicity bump luminosity star with a mass M ^ 0.8 Mq. The quantities (wrth) 
and {Dt^I K) in Fig. [3] have been scaled appropriately to reveal both our guessed functional 
dependence ffTU]) and instability condition f[T31) . At a fixed value of A"^ our solutions for the 
ratio d/Hp varying from 0.01 down to 0.0001 coincide. From Fig. [3], we conjecture that 

^mix T-th \dj V K l^/ll/ 

i.e. the thermohaline convective instability may develop {uj > 0) only if < K\N^\/Np. 
The latter condition is equivalent (ignoring factors of order unity) to (fT3l) for < K. 
If -Dh ^ K then we can neglect the term in the parentheses. In this case we obtain an 
expression for -Dmix similar to Ulrich's original formula. 

Note that an equation similar to f|T8|) can be derived directly from the dispersion relation 
f|T7|) in the same way that Ulrich used to estimate -Dmix in ffTOj) . Following him, we ignore 
the viscosity and take advantage of the fact that the thermohaline mode is the slowest one 
(cuTth ~ Tth/Tmix ^1). Therefore, we can neglect the quadratic and cubic terms in f[T7|) as 
well as terms with t:^^. Taking into account that A^^ ^ l^^l^ and A^^ ^ (TthT^)"^? we finally 
obtain 



which is reduced to ( IT8l) if z/moi < -Dh < K. 



1 - 



T 



{K + D^) \Nl\ 



(19) 



4. Observational Constraints on the Extra-Mixing Rate 



The ^He mechanism has a different underlying origin than rotational mixing, and it is 
therefore useful to re-evaluate the global implications of its operation. This is especially true 
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because we have identified two stable branches with very different timescales. We begin by 
establishing that a mixing process which occurs over too short of a timescale would have a 
dramatic impact on the thermal structure which contradicts the observational data. We also 
demonstrate that diffusion coefficients derived from empirical mixing estimates are consistent 
with the data. 

We then critically examine whether thermohaline mixing is capable of reproducing the 
data, and the answer depends critically on the assumed geometry of the fluid elements (and 
the rate at which they can achieve thermal balance with their surroundings). However, 
we can also extend the same underlying mechanism to predict luminosity trends, behavior 
in other evolutionary states, and implications for interacting binaries. Definite predictions 
emerge, and we can identify both existing conflicts and topics which will require further 
calculation. In particular, we contend that the natural expectation would be for mixing that 
is shallower in temperature and weaker for bright giants than current data indicates. The 
predictions for other types of stars are fundamentally different, and in our view the latter 
test will ultimately prove to be decisive. 



4.1. Rapid Mixing in the Adiabatic Limit 



The observed changes in the surface abundan ces of Li and C and in the ^^C/^^C ratio 



as a function of lurn inosity in LM-RGB stars (e.g.. ICharbonnel. Brown, fc Wallersteinlll998 



Gratton et aPboOoh can be used to constrain the depth and rate of extra mixing in them. If 
A log T is the difference between the logarithms of temperature at the base of the H-burning 
sh ell and at the maximum dep t h of e xtra mixing and -Dmix is the diffusion coefficient then. 



as iDenissenkov fc VandenBerg] (120031 ) have demonstrated, extra mixing in LM-RGB stars 
can be parameterized by any pair of correlated values within the close limits specified by 
AlogT ^ 0.19 and D^ix ~ 4 x 10^cm2s-\ to AlogT ^ 0.22 and D^ix ~ 8 x lO^cm^s'^ 
However, that parameterization did not take into consideration the thermal response of the 
radiative zone to mixing. 

In this paper, we use a more consistent parametric prescription by letting -Dmix be equal 
to a fixed fraction of K and allowing the temperature gradient in the radiative zone to be 
modified as prescribed by the mixing length theory relation (jl]). It is inspired by a similarity 
between equatio n (lH) and the func t ional dependence of -Dmix on K obtained for rotational 
shear mixing by lMaeder fc Meynetl (Il996l ). 



W e derive F = 0.01 for the mixing depth A log T 
( 120031 ) ■ which results in the diffusion coefficient Z^^ 



0.19 constrained by IDenissenkov &: VandenBerg 
= 0.02 7^ (assuming that I ~ d and 
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Dh K in equation [T]) . These models reproduce quite well the chemistry of LM-RGB stars 
above the bump luminosity (dashed curves in Fig. Hj). They do not noticeably change our 
models' photometric behavior and evolutionary time scale near the bump luminosity com- 
pared to models without extra mixing. This is important because photometric observations 
of LM-RGB stars, in particular the absolute V-band magnitude of the luminosity bump and 
the excess n umber of stars in the bum p, agree with predictions of standard stellar evolution 



theory (e.g.. lBjork &: Chaboyeiil2006l ). Our test computations have shown that the ratio of 
the gradients in equations (lH]-[in]) does not change much with radius and that it is roughly 
proportional to the abundance of ^He left in the mixing zone. 

On the other hand, already at F = 0.4 the bump luminosity zigzag gets so extended to- 
ward a lower luminosity (dashed curve in upper panel in Fig. [T]) , and the model star spends so 
long time following it that this peculiar behavior would sure have been noticed in photomet- 
ric studies, such as the counting of the number densities of stars as a function of luminosity 
on the RGBs of globular clusters. To be more specific, it takes about twice as long for the 
model star to make the extended zigzag as compared to the standard evolution. It should 
also be noted that the model of such a rapidly mixed star spends most of this time residing 
near the bottom of the zigzag, about 0.3 magnitude below the standard bump luminosity. 
Values of F > 0.4 would result in even more drastic changes . This behavior is similar to 



that o f models of rapidly rotating RS CVn binaries found by iDenissenkov. Chaboyer. &: Li 



( 120061 ) except that in the latter case the extended zigzag was produced by an increase of V 



caused by the stars' rotational deformations. 

The estimate of the turbulent velocity f ^ ~ gHp{A^/ fi) for the ^He-driven mixing used 
by EDL06 can be obtained from equation ([2]) if we put into it V — V = = 0. Indeed, in 
this case we can approximately consider that ~ Hp \N'^\ = gHplV^l = g H p{d In ^/dr) ~ 
gHp{Afi/ fi), at least in the vicinity of the ^He-burning shell where the mean molecular 
weight height scale is of order Hp (bottom panel in Fig. [T]). Note that in the mixing length 
theory the approximations V' = V and = are correct only in the limits of F = cxd or 
F = 0, and F^ = oo (eqs. [H-SD. 

Neglecting the influence of extra mixing on the radiative zone's thermal stratification, 
we have computed the evolution of our model also with the following diffusion coefficient: 

^mix = i:Hpv, where v'^ = ^gHp\V^\. (20) 
o o 



The factor | comes from the mixing length theory (IWeiss et al.ll2004j . Ch. 14). The depth of 
this "^He-driven" extra mixing has been determined by locating the minimum of fi above the 
H-burning shell. Outside of this point, fi increases with radius due to the ^He burning and 
mixing. We think that our prescription is in line with that EDL06 had in mind. It should 
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be noted that in our computations values of -Dmix were determined at each time step using 
a current distribution of fi that was constantly modified by extra mixing. Characteristic 
values of -Dmix obtained in this self- regulating way were of order lO^^cm^s"^. Such fast 
extra mixing is known to produc e large amounts of ^Li via th e Cameron- Fowler mechanism 



Denissenkov fc WeissI l200d : iDenissenkov fc Herwigi 120041 ) . This disagrees with the low 



(e.g. 

(often undetectable) Li abundances in the majority of LM-RGB stars located above the bump 
luminosity (compare the solid curve above log L/Lq ^ 1.8 with the observational data points 
in top panel in Fig. Hj). In accordance with EDL06, we did find a modest decrease in the 
i2q^i3q j-a^io. However, it is obvious that the observed evolutionary changes of the surface 
chemical composition of LM-RGB stars require a slightly deeper (in order to reproduce the C 
depletion) and much slower (in order to keep the Li abundance low) extra mixing than that 
advocated by EDL06. Besides, extra mixing with diffusion coefficients -Dmix ^ -ft' ~ 10®- 
10^° cm^ s~^ would bring the radiative zone to the quasi-adiabatic state (unless the fluid 
elements have flnger-like structures, which has not been reported by EDL06), which would 
cause the star to make a prolonged excursion below the bump luminosity in contradiction 
with observations. Therefore, we believe that the rapid mode originally invoked by EDL06 
does not operate, but that mild mixing is indicated. 



4.2. Mild Thermohaline Mixing in the Radiative Limit 

Because the ^He-driven thermohaline convection is expected to work in the radiative 
limit it is interesting to test if it is fast and deep enough to explain extra mixing in LM- 
RGB stars. In Fig. [5l we illustrate characteristic diffusion coefficients. In order to produce 
mixing, a physical mechanism must operate over a timescale shorter than the inflow rate. 
The dot-dashed curve shows our empirically constrained diffusion coefficient -Dmix = 0.02 K. 
For comparison, the dashed curve shows a minimum threshold diffusion coefficient -Dinflow = 
\r\Hp, where |f| is a mass inflow rate of H-rich material that flows from the bottom of 
convective envelope toward the H-burning shell (resembling a spherical accretion). The 
bottom solid curve shows a proflle of the diffusion coefficient in our unmixed bump 
luminosity model. Once mixing ensues, the //-gradient inversion spreads out over the entire 
radiative zone above the ^He-burning shell. The fin al state is illustrated with the top solid 



l ine. W e conclude that, in the prescription given by iKippenhahn. Ruschenplatt. fc Thomas 



(Il980l ) , the thermohaline convection could marginally commence and mix a narrow region in 
the vicinity of the ^He-burning shell (at r ^ 0.08 -Rq in Fig. \5^. However, the Kippenhahn 
et al. diffusion coefficient is only marginally large enough to trigger this process and is two 
orders of magnitude below the empirical value. 
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If we adopt the Ulrich prescription ffTOj) with I ^ lOd, our diffusion coefficient is raised 
by a factor of 10^ a n d the mechanism may be viable. This is the approach advocated by 
Charbonnel &: Zahru (120071 ) (incidentally, their paper was posted on astro-ph on the same 
day when we submitted the first version of our paper). 



Unfortunately, I Charbonnel fc Zahru (120071 ) do not explain how they have chosen the 
depth of mixing. We find that the depth corresponding to a minimum on the /i-profile 
(solid vertical line segments in Fig. is not sufficient to produce the observed C depletion 
(top second panel in Fig. Hj). An overshooting on a length scale of order Hp could do it 
(dotted vertical line segments in Fig. [6] are placed at a distance Hp below /imin) but then 
thermohaline "fingers" would have to penetrate down a region of higher n where they should 
experience a strong breaking. It should also be noted that the penetration of a region with 
the negative dfi/dr below /imin would reduce the average mixing rate by decreasing the slope 
of the positive d^/dr in the mixing zone. Given these uncertainties that cannot be resolved 
from first principles but instead require empirical calibrations and/or higher resolution 3D 
hydrodynamic simulations, we postpone the use of equation (fTOj) to our future paper. 



5. Concluding Remarks 



In this paper, we have shown that the ^He burning in the radiative zone of an LM-RGB 
star may drive convective fluid motions provided that their heat transport efficiency is either 
very high (the ad iabatic limit) or extr e mely low (the radiative limit). Confirming the con- 
clusions made by ICharbonnel fc ZahnI (120071 ) , we identify the mixing in the radiative limit 
with thermohaline convection and we note that this convection would have a sufficiently 
high rate to explain the observed mixing pattern in LM-RGB stars only if fluid elements 
could travel over length scales exceeding their diameters by a factor of 10 or more. However, 
we also find that thermohaline convection may be suppressed by horizontal turbulence if its 
associated diffusivity -Dh ^ 3fr|V^|/(Vad — Vran)- Such va l ues o f for rot at ion- induced 
horizontal turbulence have bee n obtained by iPalacios et al.l (I20061) who used a prescription 
for estimating proposed by lMathis. Palacios. fc ZahnI (120041 ). Although this prescription 
may be considered as a very primitive approximation to a complex physical phenomenon, a 
similar heuristic approach has been used to successfully m odel mixing and angular momen 



tum t ransport in ra diative zones of massive MS stars (e.g. jTalon fc Zahrull997l : iTalon et al. 



19971 : iMaederl |2003| ). Of course, we recognize that rigorous 3D hydrodynamic simulations 
have yet to be done to support or refute these heuristic models. In addition to the theoretical 
issues above, there are significant empirical challenges for an explanation that relies solely on 
thermohaline convection. The same physics should consistently be applied to other phases 
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of evolution or situations where n inversions occur. 



If extra mixing in RGB stars is really driven by ^He burning then it should die out by the 
end of the RGB evolution because of the ^He exhaustion. In this case, the ^He-driven extra 
mixing could not resume working in the same stars on the asymptotic giant branch (AGB). 
So, we would expect the absence of observational signatures of extra mixing in low-mass 
(M < 2 Mq) AGB stars unless the mixing in them is of a different nature. However, given 
the similarities in their depth and in the structure of radiative zone where they operate, it is 
unlikely that the RGB and AGB mixing have different physical mechanisms. Contrary to this 
prediction, there are observation al data indicating the presence or necessity of operation o f 



extra mixing in these stars (e.g.. iNoUett. Busso. fc Wasserburg||2003l : iMasseron et al.ll2006l ). 



Moreover, to comply with observations, the AGB mixing has to penetrate close enough to 
the H-burning shell to dredge up material processed in the CN-cycle, like in RGB stars, 
mimicking the convective hot-bottom burning that occurs in more massive AGB stars. 



Stancliffe et al.l (120071 ) have used the Kippenhahn et al. prescription (Q to model ther- 



mohaline mixing in a metal-poor low-mass MS star accreting wind material from its AGB 
binary companion enriched in He and C. Such accretion is believed to be the primary process 
responsible for the formation of the so-called carbon-enhanced metal-poor (CEMP) stars. 
They have found that thermohaline convection mixes almost 90% of the star within about 
10^ years after the accretion. On the other hand, the RGB mixing pattern can be repro- 
duced only if the diffusion coeffic ient given by equation is increased by a factor of 10^ 
to 10^ (jCharbonnel fc Zahij 120071 : also see our Fig. In this case, thermohaline mixing 
in CEMP MS stars would dilute the accreted material on a much shorter time scale of or- 
der 10^-10^ years. Unless most of the CEMP stars accreted substantial fractions of their 
initial masses, their r ather high frequency /cemp ^20% among very metal-poor stars (e.g., 
Lucatello et al.l 2006 ) would look surprising. Furthermore, both Lucatello et al. ( 2006 ) and 
Aoki et al.l (120071 ) have found anti-correlations between [C/H] (and [(C+N)/H]) and lumi- 
nosity for Ba-enhanced CEMP stras spanning over three orders of magnitude in L that they 
interpreted as an evidence of dilution of the envelope material in the accreting companion. 
Thermohaline convection on a time scale of order 10^ - 10^ years would mix the CEMP MS 
stars almost instantaneously and well before their luminosity would begin to increase due 
to the core H exhaustion. In that case, the mentioned anti-correlations could not have ap- 
peared. Also note that even on the lower RGEl Ba-e nhanced CEMP stars have quite low 
carbon isotopic ratios (^^C/^^C < 20; iRyan et al.ll2005l ) in a striking contrast with the values 
of ^^C/^^ C^ lOOQ tha t the low-mass AGB stars are predicted to return to the interstellar 
medium (jHerwig]|2004l ). Extra mixing (in the low-mas AGB stars) could easily resolve this 
discrepancy. 
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The problem of mixing i n CEM P MS and RGB stars has recently been addressed by 
Denissenkov &: Pinsonneaultl (l2007l ). Particularly, they have shown that the first dredge- 
up dilution of CN enrichment in CEMP stars relative to their MS precursors is indeed a 
plausible explanation of the observed anticorrelation of [N/Fe] with log L/Lq and that it 
contradicts models that rely on efficient thermohaline mixing indu ced by sma l l /x gr adients 
in red giants. This result has independently been confirmed by lAoki et al.l (120081 ). The 
suppression of thermohaline convection by rotationally driven horizontal turbulence may 
explain its reduced efficiency in MS CEMP stars. 

Another potentially serious problem for the ^He-driven thermohaline convection could 
be to explain available observational evidence of enhanced extra mixing in LM-RGB stars. 
Firstly, observations show that in some globular clusters the anti-correlated abundance vari- 
ations of C and N in red giants become larger when stars approach the RGB tip. More- 
over, extremely large values of the N abundance in some of these stars indicate the dredge- 
up of material in w hich not only C but also a fraction of O has been converted into N 
(jSmith et al.ll2005al ). Secondly, at least in the globular cluster M13, the rela tive number of 



uppe r RGB stars with the 0-Na anti-correlation increas e s with luminosity (jJohnson et al. 



20051 ). Thirdly, iDenissenkov. Pinsonneault fc Terndrupl (120061 ) have shown t hat the ^^F 



abundance variations found in bright red giants of the globular cluster M4 by ISmith et al. 



(l2005bl ) may also require that extra mixing in them to operate much faster and somewhat 
deeper than in LM-RGB stars immediately above the bump luminosity. It is difficult to in- 
terpret these data by the ^He-driven mixing because its efficiency should decline toward the 
RGB tip in proportion as '^He gets depleted. Unfortunately, the question of evolutionary Na, 
O, and ^^F abundance variations in globular-clusters RGB stars is still a matter of debate 
from the observational point of view. Therefore, we consider them as a potential rather than 
a real problem for the ^He thermohaline convection. 

A similar problem is encountered when one tries to understand the phenomenon of Li- 
rich giants. There are convincing arguments that high Li abundances in these LM-RGB 
stars are produced via the Cameron- Fo wler mechanism that also requires enhanced extra 
mixin g with -Dmix ~ 10^° -10^^ cm^ s~^ (IDenissenkov &: WeissI l2000l : iDenissenkov fc Herwig 
20041 ). It should be noted that most of the Li-rich giants are located above the bump lumi- 
nosity (jCharbonnel &: Balachandranll2000l ). Besides, their proportion among rapid rotators 
(f sin i > 8kms~^) is ~50% which is considerably larger than ~2% of Li- rich stars amori g 
the much more common slowly rotating (v sin i < Ikms"^) K-giants (IDrake et all 120021 1. 
It is no t clear why the ^He-driven mixing would be enhanced in fast rotators. Oppositely, 
Canutd (119991 ) argues that a larger shear due to differential rotation should decrease the 
efficiency of mixing by thermohaline convection. A larger shear would also intensify horizon- 
tal turbulence, thus hindering thermohaline convection even stronger. So, the Li-rich giants 
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seem to support the hypothesis of rotation-induced mixing rather than that of thermohahne 
mixing. This may not necessarily be rotational shear mixing that has already been criti- 
cized. Instead, rotation may deposit its kin etic energy to mix ing less directly, e.g. through 



generation of buoyant magnetic flux tubes (IBusso et al.l 120071 ). 



For the adiabatic limit, the predicted evolutionary changes of the surface composition of 
LM-RGB stars disagree with observational data. Besides, high F's would bring the radiative 
zone to the quasi-adiabatic state which would result in a photometric behavior of the RGB 
star inconsistent with the observed one. 

In the vicinity of /Xmin, our empirically constrained diffusion coefficient has values of order 
-Dmix ~ 10^-10^ cm^ s~^ (dot-dashed curve in Fig. [5]). If we assume that in real LM-RGB 
stars extra mixing is produced by thermohaline convection whose fluid elements have a ratio 
of ~ 10, where I does not exceed the local pressure scale height Hp ~ 0.02 i?©, then we 
can estimate the elements' characteristic velocities v ~ 3Dmix/^ ~ 2 x 10~^-2 x 10~^ {Hp /I) 
cms~^. EDL06 have found velocities of order 5 x 10^ cms~^ in their 3D red giant model. 
Those would correspond to fluid elements with diameters from 6 to 60 cm! Interestingly, such 
small fluid elements would actually be optically thin because the photon mean free path is 
~ 1 cm in this environment. However, we do not believe that EDL06 could resolve such small 
finger-like structures. As we mentioned, they have reported fiuid element displacements of 
order / ~ 10^ cm. Given the discussed inconsistencies in modeling extra mixing in LM- 
RGB stars with the ^He-driven convection, we conclude that a different mechanism is worth 
searching for. It is also obvious that higher resolution 3D hydrodynamic simulations of 
the ^He-driven mixing are needed to understand what EDL06 have actually witnessed. In 
particular, we qoute here an issue raised by the referee, who seems to be an expert in the 
field. "The EDL06 calculations had no mechanism to simulate the turbulent cascade on scales 
smaller than the zoning (a subgrid scale model, such as an eddy viscosity treatment), and 
yet the calculations did not numerically blow up. To me this says that the finite difference 
expressions of the Djehuty code are themselves quite diffusive. Thus, I am surprised that this 
incredibly small mean molecular weight inversion could generate significant motion without 
being squelched by the numerical diffusion. My conclusion is that there are some serious 
issues that must be addressed about the numerical behavior of the 3D calculations." 



We thank the anonymous referee for useful comments and suggestions that helped us to 
improve the manuscript. We acknowledge support from the NASA grant NNG05 GG20G. 
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Fig. 1. — Upper panel: evolutionary tracks of a 0.83 M0 model star (the initial H and He 
mass fractions are X = 0.758 and Y = 0.24) near the bump luminosity (solid curve — for 
r = 0.01; dashed curve — for F = 0.4). Cross marks the end of the first dredge-up. Bottom 
panel: profiles of the mean molecular weight in the radiative zones of our unmixed red giant 
models locations of which are shown with the same symbols in the upper panel. Depressions 
of fj, are caused by the reaction ^He (^He, 2p)^He. 




Fig. 2. — Profiles of tlie ratio of tfic thermal and molecular difTusivity (solid curve) and of 
the quantity Kq = cm^s~^ (dashed curve) in the radiative zone of our 0.83 Mq bump 

luminosity model. 




Fig. 3. — Solutions of the dispersion relation f|T7|) as a function of the horizontal turbulent 
diffusion coefficient. Both the only real root (curth) and the parameter [Dh/K) have been 
scaled appropriately to reveal the dependence ffTOl) and the instability condition f|T3|) . Plotted 
are the solutions for the fluid element diameter d = 0.01 Hp and three values of A^"^ = —10"'' 
(this is a characteristic value for ^He burning in low-metallicity bump luminosity stars with 
M pa 0.8 Mq] solid line), —10"® (short-dashed line), and —10"^ (long-dashed line). 
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Fig. 4. — Comparison of the observational data from iGratton et al.l ( 20001 ) for field metal- 
poor (—2 < [Fe/H] < —1) low-mass stars (circles) with results of our computations of the 
evolution of the 0.83 Mq star with the ^He-driven mixing (solid curves, equations [20]) and 
extra mixing with the rate -Dmix = 0.02 K and depth A logT = 0.19 (dashed curves). Dotted 
lines show predictions of the standard theory. For further details, see text. 
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Fig. 5. — Dot-dashed curve - the empirical profile of Z^mix = 0.02 7^; dashed curve - a 
minimum threshold profile -Dinflow = \r\Hp that any Z^mix must exceed; bottom solid curve 
- the -Dxipp (thermohaline convection, equation [9]) profile in our unmixed bump luminosity 
model (a hump ai r ^ 0.08 -Rsun is produced by a local increase of |V^| in the ^He-burning 
shell); top solid curve - -Dxipp in a model in which mixing with -Dmix = 0.02 K has spread 
out the /i-gradient inversion over the entire radiative zone. 
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Fig. 6. — Profiles of the mass fractions of ^He and CNO elements in the vicinity of the 
H-burning shell in 0.83 Mq bump luminosity models with the heavy-element mass fractions 
Z — 0.002, and Z — 0.0001. Vertical solid line segments show locations of the minimum on 
the /x-profile. Dotted line segments arc placed one pressure scale height below ^min- Without 
overshooting, the depth of the ^He-driven thermohaline convection would be at the locations 
of the solid segments. There would be no evolutionary C depletion in this case, contrary to 
observations. 



